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A Possible Mechanism of the Pseudogap in Organic Superconductor 
Based on the Superconducting Fluctuation 

Takanobu JuJoQ and Kosaku Yamada 

Department of Physics, Kyoto University, Kyoto 606-8502 
(Received March 19, 1999 ) 

The explanation of the anomalous behavior in K-type (BEDT-TTF)2X which was revealed by 
the nuclear magnetic resonance experiments is presented. We calculate the electronic properties 
by using the one-loop approximation for the superconducting fluctuation. The decrease of the 
density of states due to the large damping effect is obtained and the one-particle spectrum has 
a pseudogap structure around the Fermi level. It is found that these behaviors are peculiar 
to the quasi two-dimensional system and are enhanced in the incoherent metal. The similarity 
between the cuprates and this compounds is discussed and an experiment to check this proposal 
is suggested. 

KEYWORDS: organic superconductor, superconducting fluctuation, pseudogap, one-particle spectrum, quasi two- 
dimensional system, incoherent metal 
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Organic compounds arc known to show various phases 
including superccmductivity under controlled tempera- 
ture and pressurcEP In these compounds K-type (BEDT- 
TTF)2X (abbreviated as k-(ET)2X hereafter, and 
X=Cu(NCS)2, Cu[N(CN)2]Br, Cu[N(CN)2]Cl) has at- 
tracted much attention because it has the highest su- 
perconducting transition temperature. The characteris- 
tic properties of k-(ET)2X are that the two-dimensional 
tight-binding approximation well explains the Shubnikov- 
de Haas experimentso-* and the superconducting phase 
neighborSjthe antiferromagnetic phase under a controlled 
pressure.™ About the properties of the superconductivijti^ 
of K-(ET)2X, the nuclear magnetic resonance (NMRp-' 
and the specific heatla^ experiments of «;-(ET)2X suggest 
that the superconducting gap has line nodes. Above the 
superconducting transition temperature (Tc), k-(ET)2X 
shows the anomalous pppperties, which were revealed by 
the NMR experiments.B' The striking character of the 
anomalous properties is that the uniform magnetic sus- 
ceptibility and l/TiT (Ti is the spin lattice relaxation 
rate and T is temperature) are not independent of tem- 
perature dependence but decrease with decreasing tem- 
perature below T ~ 50K. (This phenomenon may be 
called the spin gap or the pseudogap.) 

Theoretically, the importance of the electron corre- 
lation in «;-(ET)2X is suggested because its dimerized 
structure makes the system half- filled Jj-' The mechanism 
of superconductivity of k-(ET)2X is investigated using 
several approaches like the fluctuation exchange approx- 
imation lELEX),*'^^ the third-order perturbation theoiy 
(T0PT)E3) and the quantum Monte Carlo simulation.E3 
From these studies, it is found that this superconductiv- 
ity is a spin fluctuation mediated one and the symmetry 
of the Cooper pair is c?-wave in this compound. How- 
ever, for the anomalous properties above Tc, there has 
been no theoretical proposal about the origin of these 



behaviors until now, and the aim of this paper is to pro- 
pose a possible mechanism of the (pseudo-)spin gap in 
K-(ET)2X. 

/pijET)2X has a strong anisotropy in the conductiv- 
ity,EJ that is, it can be considered that this compound 
is a quasi two-dimensional system. From FLEX and 
TOPT, it is also known that Tc of k-(ET)2X is rather 
high when Tc is scaled with the bandwidth. Regard- 
ing the properties of k-(ET)2X stated above, it can be 
considered that the superconducting fluctuation plays an 
essential role in the appearance of the (pseudo-) spin gap. 

There are many theories on the pseudogap in cuprates 
on the basis of the superconducting fluctuation, but these 
calculations are restricted to a two-dimensional (2D) case^^'^^^ 
and the effect of the dimensionality has not been consid- 
ered yet. 

In this Letter, we investigate the electronic properties 
for several values of the anisotropy parameter and for 
those of the attractive force by using the one-loop ap- 
proximation for the superconducting fluctuation. From 
the nature of the one-loop calculation, our interest doesn't 
lie in the effect of the fluctuation on Tc, but in the elec- 
tronic properties above Tc. (Therefore Tc is the mean 
field transition temperature in our discussion.) It is re- 
vealed that, owing to the superconducting, fluctuation 
the damping rate of quasiparticles increases and then 
the density of states (DOS) at the Fermi level decreases; 
this possibly explains the anomalous behavior. We also 
obtained the result that the one-particle spectrum has a 
pseudogap structure. This behavior is expected to be ob- 
served in the angle-resolved photoemission spectroscopy 
(ARPES) experiments. It is found that this pseudogap 
behavior is a characteristic property of a quasi 2D sys- 
tem such as K-(ET)2X. Furthermore, this behavior is en- 
hanced in the incoherent metallic region where the spec- 
trum is thermally broadened. 

To observe the effect of the dimensionality, we consider 



* E-mail: jujoOton. scphys.kyoto-u.ac.jp 



2 



Takanobu JujO and Kosaku Yamada 



the foUowing Hamiltonian, 

^^Y.'k'i,a^k,a- E ak,k'^Jq-kf'q-k\i''k\r 
k,a k,k',q 

(1) 

Here ej^ is the noninteracting energy dispersion, 

ej^ — —2t{coskx+cosky)—2t'cos{kx+ky)—2tzcos{2kz)~ii. 

is the chemical potential, and we adopt the dimer 
model. Therefore fi is determined so as to fix the particle 
number per site to 1. We take t as the unit of energy, 
and fix the ratio t' /t to 0.70 as a realistic value. It should 
be noted that we include the inter-plane hopping term in 
the energy dispersion. The essential point to study the 
effect of the dimensionality, is not the form of this term 
but the ratio between t, t' (intra-plane hopping) and tz 
(inter-plane hopping). We take the separable d-wave in- 
teraction as the attractive force, that is, ^/ — 94'^^]^' j 
= cos{kx) - cos(fcy)). 

On the basis of this model, we take the superconduct- 
ing fluctuation into account by the one-loop approxima- 
tion, which is taking T-matrix up to the flrst order in 
calculating the self-energy. The self-energy is written as 
the following form; 



s(fc) 



-E 



Tk,k(i)'^o{q- k)- 



(3) 



Here, N is the number of sites, Go is the bare Green's 
function and j^/ (g) is the T-matrix; 



where D{q) is the boson propagator; 



iEfe0|Go(fc)Go(q-A:)- 



(4) 



(5) 



Here, k is (fc, e„) (e„ is the fermionic Matsubara fre- 
quency e„ = (2n + l)7rT, n is integer), and q is {q,uj„i) 
{uJm is the bosonic Matsubara frequency ujm = 2rmrT, 
m is integer). The numerical calculations are performed 
with 32 X 32 X 32 meshes in the first Brillouin zone, and 
with 512 Matsubara frequencies. This number of fre- 
quencies ensures reliable results down to T = 0.003. 

Now, we show the results obtained by our calculation. 
First, we see the effect of the dimensionality by calcu- 
lating the electronic properties for different values of ■ 
The calculated self-energies for the inter-plane transfer 
tz = 0.10 and 1.00 are shown in Fig. 1. The momentum 
of the self-energies is fixed to the Fermi momentum situ- 
ated on the line connecting F point (0, 0, 0) and M point 
(tTjO, 0). At this point called A (see Fig. 5(a)), the at- 
tractive force takes the strongest value due to the dx2_y2- 
wave symmetry. For the isotropic case (tz = 1.00), the 
behavior of the self-energy is the same as that of the con- 
ventional Fermi liquid. That is, — Iml](fc,a;) has a local 
minimum at w = and the slope of Ke'E{k,uj) around 

= is negative. These properties oLthe self-energy 
indicate the well-defined quasiparticle.Ej^ However, for 
the quasi 2D case {tz = 0.10), the behavior of the self- 
energy is opposite to that of the isotropic case. The peak 




Fig. 1. The self-energy for = 0.10 (Tc = 0.337) and = 1.00 
(Tc = 0.195). The values of g and temperatures are shown in the 
figure, (a) The real part of the self-energy (b) The imaginary 
part of the self-energy. 



of —lmY,{k,uj) at ti; = indicates that the quasiparticle 
is not well-defined, and the positive slope of ReS(fc,a;) 
around uj ~ indicates that the spectral weight of quasi- 
particles near the Fermi surface decreases. Therefore the 
DOS at the Fermi level is expected to be suppressed. 
The fact that there exist these anomalous behaviors of 
the self-energy behind the precursor of the_auperconduc- 
tivity was first pointed out by Janko, et a/o^ In addition 
to this fact, it is pointed out here that these anomalous 
properties are seen only in the quasi 2D system. 

To observe the relationship between the self-energy 
and the DOS, the temperature dependences of p(0) (the 
DOS at the Fermi level) and — ImE(fcF,0) are shown 
in Fig. 2. It can be seen that fo r the quasi 2D case 
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Fig. 2. The temperature dependences of the DOS at the Fermi 
level p{0), and the absolute value of the imaginary part of the 
self-energy at the FS and the Fermi level — ImS(fcp,0). Here fcp 
is at the point A. "ise" in the figure means the imaginary part 
of the self-energy. Tc is the same as in Fig. 1. 



{tz = 0.10) p{Q) begins to decrease far above Tc with the 
increase of the damping rate of quasiparticles near the 
Fermi level, — ImI](fcF,0). On the other hand, for the 
isotropic case {tz = 1.00) p(0) increases down to Tc be- 
cause the system acquires coherency as the temperature 
is decreased, which is supported by the decrease of the 
damping rate of quasiparticles, — ImE(fcFjO). Although, 
as discussed below, two other factors are required for 
the quantitative comparison with experimental results, 
the decrease of p(0) over a wide temperature range sug- 
gests that the precursor of the superconductivity can be a 
candidate mechanism which gives rise to the anomalous 
behavior in k-(ET)2X, because the quantities showing 
the anomalous behavior in the experiments are related 
to the DOS at the Fermi level. 

The above discussion is described analytically by as- 
suming that the boson propagator is written in the fol- 
lowing form; 
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D{q,Lu) 



1 



(6) 



M - ia,uj + ■ 

Here, M is the mass term which describes how far the 
system is from the superconducting state and is written 

as M = i - Efc 0| ^ Af = at T = T,. 

With the use of this form for the boson propagator, 
— ImE(fcF, 0) can be written in the foUowing form; 



- ImI](A;F,0) 
for a 3D system, and 



- ImI](fcF,0) 



log(H-%L), (7) 



kF 



M 



C9c 



,tan-i(^) 



(8) 



2ttvf£,vM VM 

for a 2D system. Here, Qc is the cutoff momentum of the 
boson, and v-p is the vefocity of the electron at the Fermi 
surface (FS). 

From these equations, it can be seen that for a 2D 
system — ImE(fcF,0) is large over a wide temperature 
region above Tc. The interpolating formula between 3D 
and 2D cases can be written following the discussion by 
Takahashi,E3' 



found in k-(ET)2X which is a quasi 2D system. On the 
other hand, for large tz, Tc is low because of the large 
bandwidth. Therefore, we studied the spectrum for large 
g in the case of large tz- However, we saw no gap-like 
behavior in this case either. Therefore we can conclude 
that the gap-like behavior is a characteristic property of 
the quasi 2D system. 

So far, we have fixed the magnitude of the attractive 
force g. Next, we will see whether the pseudogap behav- 
ior can be observed for all values of g. The calculations 
are performed by determining the value of M for which 
the pseudogap opens for various values of g. Using this 
method, we can calculate the one-particle spectrum at 
the same point away from the superconducting state for 
various values of g. The results ar e shown in Fig. 4. It 
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Fig. 4. The relationship between the attractive force (g) and the 
temperature (Tpg) at which the pseudogap opens in the one- 
particle spectrum at point A. The value of tz is shown in the 
figure. M is the mass-term. 



-ImE(fcF,0) 



fcp 



fcF 



167ruF^^e 



log(l 



M 



[tan-i(^)-tan-i(^)]. (9) 



M 



Here, e tz/t is the parameter which characterizes the 
degree of the anisotropy. 

We also calculated the one-particle spectrum {A{k, uj) = 
—lmG{k,uj)/TT), and it is found that for the small inter- 
plane transfer tz the spectrum has a pseudogap structure 
around uj = 0. The relationship between tz and the tem- 
perature (Tpg) at which the pseudogap opens at point 
A, is shown in Fig. 3. Tpg is obtained in a manner such 
that the one-particle spectrum is calculated at intervals 
of 0.001 at a temperature which is less than one percent 
of Tc, and Tpg is the temperature at which it is found to 
be suppressed around uj = 0. It is seen that for the small 
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Fig. 3. The relationship between the anisotropy parameter (tz) 
and the temperature (Tpg) at which the pseudogap opens in the 
one-particle spectrum at the point A. The value of g is shown 
in the figure. M is the mass-term defined above. 



inter-plane transfer tz the gap-like behavior can be seen 
at rather high temperatures above Tc.EZP Therefore, the 
pseudogap in the onc-particlc spectrum is expected to be 



can be seen that the pseudogap behavior is enhanced for 
large g. This result is interesting because the self-energy 
is independent of g when the self-energy is written as a 
function of M and therefore it is expected that the gap- 
like behavior appears for the same value of M . However, 
our results show that this assumption is incorrect. The 
behavior in Fig. 4 is explained by the following argu- 
ment. The self-energy is written from eq. (8), 

- ImS(feF, 0) cx y%{^fJ=. (10) 



^Tf' 

Here Tf is the Fermi degeneracy temperature. For large 
g, Tc is high and therefore the enhancement factor {T/TyY 
is large. Therefore, in the case of high Tc, the damping 
effect is more enhanced even for rather large values of M 
than in the case of low Tc. In other words the thermally 
broadened spectrum is more sensitive to the supercon- 
ducting fluctuation. 

This situation is considered to be realized because the 
superconductivity of k-(ET)2X appears only near the 
Mott-insulating phase and it is likely that the Tf is very 
low owing to the strong electron correlation. From the 
above discussion, we propose that point A and its symmetry- 
related points of the Fermi surface must be destructed in 
this pseudogap region due to the strong fluctuation of the 
superconductivity. This will be checked by the ARPES 
experiments. The FS and the one-particle spectrum at 
point A are shown in Fig. 5. In Fig. 5(a), the FS of the 
single-band in our model is shown. The equivalence of 
this single band model and the two-band model (the hop- 
ping parameter t takes two kinds of values) is discussed 
in ref. 8. From Fig. 5(b), it can be seen that the one- 
particle spectrum has a pseudogap structure above Tc. 
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Fig. 5. (a) The Fermi surface for = 0.00. The meaning of 
point A is the same as above, (b) The temperature dependence 
of the one-particle spectrum. The values of g and are shown 
in the figure. In this case Tc is 0.339. 



This behavior disappears at the points of the FS away 
from point A and its symmetry-related points because of 
the d2.2„j,2-wave symmetry of the attractive force. 

Finally, we comment on the two factors which are not 
included in our discussion but which are required to be 
considered for an accurate quantitative description. One, 
is that we do not take into account the effect of the fluctu- 
ation in the determination of Tc, i.e. the self-consistency 
requirement. For a quasi 2D system such as k-(ET)2X, 
the fluctuation plays an essential role in determining T^. 
However, our conclusion that the pseudogap is peculiar 
to the quasi 2D system is not affected because our ar- 
gument is basing on the mass term M which measures 
the extent to which the system is away from the super- 
conducting state. The other factor is the effect of the 
Coulomb repulsion. This is also important for obtaining 
accurate quantitative results in the system which neigh- 
bors the Mott insulating phase. However, if the adia- 
batic continuity is realized, the effect of the Coulomb 
repulsion is only the rescaling of various quantities, such 
as the Fermi degeneracy temperature. Justification of 
the adiabatic continuity in the metallic phase near the 
Mott-insulating phase is an important problemjl3^ and 
will be discussed elsewhere. 

In summary, we have calculated the electronic proper- 
ties of the quasi 2D model by considering the supercon- 
ducting fluctuation. The relationship between the DOS 
and the damping rate is shown and it is found that the 
decrease of the DOS is caused by the large damping rate 
which is enhanced by the superconducting fluctuation 
in the quasi 2D case. We have also calculated the one- 
particle spectrum, and the pseudogap behavior over a 
wide temperature range above Tc is obtained only for the 
quasi 2D case. It is also found that this behavior is char- 
acteristic of a rather incoherent metal which has a broad 
spectrum around the Fermi level and is sensitive to the 
fluctuation. These conditions for the appearance of the 
pseudogap can be applied to k-(ET)2X. Therefore, the 
anomalous properties found by the NMR experiments are 
considered to be caused by the superconducting fluctua- 
tion. The strong two-dimensionality and the incoherent 
metallic properties are also the properties of cuprates in 
the under-doped region, so our argument can be applied 
to this material and supports the expectation that the 
anomalous behaviors of the cuprates and k-(ET)2X have 
the same origin.®' ARPES experiments are desired not 
only to check the validity of our proposal, but also to un- 
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derstand the property of pseudogap in k-(ET)2X more 
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